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Improving the locality of the overlap Dirac operator via approximate 
solutions of the Ginsparg- Wilson relation 

David H. Adams'' * 

^Institimt-Lorentz for Theoretical Physics, Leiden University, 
Niels Bohrweg 2, NL-2333 CA Leiden, The Netherlands 

We determine a free field hypercubic lattice Dirac operator which is optimally close to satisfying the Ginsparg- 
Wilson relation. Inserting this operator into the overlap formula, we show that the analytic locality bound on the 
resulting overlap Dirac operator is substantially stronger than in the standard case. The improvement generally 
persists in gauge backgrounds when the plaquette variables are all close to unity. 



The overlap Dirac operator provides a formula- 
tion of lattice QCD with exactly massless quarks 
Actually there are many overlap operators: 
for each choice of ultra-local lattice Dirac opera- 
tor D an overlap Dirac operator D^y is obtained 
by inserting D into the overlap formula: 



D,., = -{l + A{A*Ay 
a 



-1/2N 



A^D-^ 



(1) 



Here a is the lattice spacing and m is a param- 
eter which controls topological properties of Doy 
and the number of fcrmion species described by 
the corresponding lattice fermion action. In the 
standard overlap operator, D is taken to be the 
Wilson-Dirac operator D^j . 

The locality issue for Dov is nontrivial due to 
the inverse square root which 
also causes difficulties for numerical implementa- 
tion. To ensure that the lattice theory is in the 
right universality class to reproduce continuum 
physics, Dov should be exponentially-local. An 
exponential-locality bound was derived in |2]: 



\\{A*A)-^/^{x,y)\\ < Ke-'^^^-y^^^ 



(2) 



where 



-y\ = \x^^-y^^\ 

imum distance in units of lattice spacing between 
lattice sites coupled by D as measured by this 
norm (e.g. I — 1 for the Wilson-Dirac operator). 
Then the decay constant in (0) is 



Set / to be the max- 



0/21 



(3) 
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with 6 determined by the condition number C: 



= log ( 



^max{A* A) 



(4) 



^max i^min) ^rc thc maximum (minimum) eigen- 
values. Moreover K = \min{A*A)~^/-^ (this and 
(j2Jl-l0J were not written out explicitly in |2] but 
follow easily from the formulae derived there). 
To obtain a locality bound independent of the 
gauge field, gauge field-independent bounds < 
u < A* A < V < oo are required; then Xmin and 
^max can be replaced by u and v, respectively, 
in the preceding. An upper bound v can easily 
be derived via triangle inequalities. However, to 
get a nonzero lower bound u some restriction on 
the lattice gauge fields is generally required, since 
gauge backgrounds in which A has zero-modes do 
exist in general (reflecting topological properties 
of Dov)- In the standard Wilson-Dirac case this 
can be done by imposing an admissibility condi- 
tion I2I3I4I 



\1-U{p)\ \ < e V plaquette p. 



(5) 



The currently sharpest lower bound is 1 — 6(2-1- 
V2)e < A'^Avj (in 4 dimensions and with m ~ 1) 

m 

In |S] the possibility of improving locality and 
convergence properties of the overlap Dirac oper- 
ator by taking the input D to be an approximate 
solution of the Ginsparg- Wilson (GW) relation 
ini was pointed out. This is based on the ob- 
servation that if D is an exact GW solution, i.e. 
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75 -D + Z575 = ^Djc^D, then the overlap formula 
simply gives this operator back again: = D. 
(To see this, note that the GW relation is equiv- 
alent to A* A = where A = D - It 
is known that ultra-local operators cannot ex- 
actly satisfy the GW relation |7]; but approx- 
imate solutions are possible, and for these the 
preceding observation implies Dov ~ D, indicat- 
ing that Dov is "close to being ultra-local" . This 
heuristic reasoning needs to be treated with cau- 
tion though, because if D is very close to satis- 
fying the GW relation then it must be "almost 
non-ultralocal" . In this note we attempt to get 
a more precise analytic understanding of the sit- 
uation by considering the improvement that can 
be achieved in the locality bound If I? is 
an approximate GW solution, i.e. A* A « (^)^j 
then Xmin ~ ^max ~ (■^)^; hcncc the condition 
number C is close to unity and consequently the 
parameter in is very large. However, the ef- 
fect of the the increase in in the decay constant 
(O must be at least partly offset by an increase 
in I , since the closer D is to satisfying the GW re- 
lation the longer its range must be. Therefore, to 
improve the locality bound on the overlap Dirac 
operator, we are led to look for short range ap- 
proximate solutions to the GW relation to use as 
input in the overlap formula. 

A natural arena in which to look for such oper- 
ators is the class of hypercubic operators. These 
are the lattice Dirac operators which couple sites 
within the same lattice hypercube. For such oper- 
ators the parameter I in Q equals the spacetime 
dimension d, so in 4 dimensions an increase in 9 
by more than a factor of 4 is required in order 
to achieve a locality bound that is stronger than 
the one for the standard overlap operator. We 
begin by considering the free field case. Define 
the hermitian operators S"^ = 51 (7+^ — T^) and 
C^, = ^{T+^+T_^) where (T_^) are the usual 
forward (backward) parallel transport operators. 
Then a general free field hypercubic operator can 
be written in terms of the free field 5'^'s and Ci,'s 
as IE] 



and A = Ep=o2^\E 



^ = -(7>M + A) 



(6) 



• • • Cvr, 



ui< - <vp - ■ ■ Here 
\d are coupling param- 



Ki, . . . , Krf and Ao 
eters; they are the most general couplings 
allowed by the lattice symmetries. The re- 
quirements of correct formal continuum limit 
{D — > 7^9^) and vanishing bare mass are equiv- 



alent to the constraints X]p=i 2^ 



d-l 
p-l 
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and Ep=o 2P p \ = «• These can be used to 
eliminate ki and Aq; then D contains 2d — 1 free 
parameters K2, ■ • • , i^d] Ai, . . . , A^. (The Wilson- 
Dirac operator, which contains just one free 
parameter -the Wilson parameter r - is obtained 
by setting m = 1/2, Aq = dr, Ai = — r/2 and 



A,, 



for p > 2.) 



A number of free field hypercubic operators 
which approximately satisfy the GW relation are 
already known and have been discussed in 
the truncation of the standard overlap operator, 
the truncated Fixed Point (FP) operator (the 
FP operator is a GW solution which is known 
explicitly only in the free field case |9j), and 
the "GW-improved" operator obtained by insert- 
ing the truncated FP operator into the over- 
lap formula and then truncating again. How- 
ever, one can attempt to do better still by con- 
sidering the condition number Q as a func- 
tion C — C(k2, ■ • ■ , Kd; Ai, . . . , Ad) on the pa- 
rameter space of free field hypercubic operators 
and finding the point (s) at which C has a mini- 
mum. We did this numerically, using the expres- 
sion C = {A*A{k)} I {A*A{k)} where 

A* A{k) is the free field momentum representa- 
tion of A* A. For concreteness we fix the space- 
time dimension to be c? = 4 and set m = 1 in l(T]l. 
i.e. a A = aD — 1. The lattice spacing a drops 
out in the expression for C. When calculating a 
minimum of C(k2, . . . , A4) numerically one must 
specify an initial point (actually two initial points 
for our Mathematica calculation). We tried var- 
ious choices of starting points corresponding to 
the truncated FP operator and other free field 
hypercubic operators considered in In each 
case the the numerical calculation converged to 
the same minimum of C{k2, ■ ■ ■ , A4) at the point 
in parameter space listed in the last column of 
Table Thus we have found a new free field 
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Trunc. FP 


Optimal 






0.03208 


0.03230 






0.01106 


0.01107 




K4 


0.00475 


0.00613 




Ai 


-0.06076 


-0.06191 




A2 


-0.03004 


-0.03027 




A3 


-0.01597 


-0.01576 


Tnhlp 1 


A4 


-0.00843 


-0.00778 


The coupling parameters of the new 


"optimal 



free field hypercubic operator. For comparison the 
parameters of the truncated FP operator are also 
listed. 







^max 


C 


e 


Wilson-Dirac 


1 


49 


49 


0.14 


Trunc. FP 


0.919 


1.023 


1.113 


0.45 


GW-improved 


0.932 


1.032 


1.107 


0.46 


Optimal 


0.938 


1.005 


1.072 


0.51 



Table 2 

^min{A*A), X„iaxiA*A), C and for the free field 
Wilson-Dirac operator, the truncated FP and its 
GW-improved version, and the new "optimal" hy- 
percubic operator. 

hypercubic operator which is optimaUy close to 
satisfying the GW relation, at least compared to 
other operators in a neighbourhood of the trun- 
cated FP operator. 

In Table |21 we list Xmin , Xmax and the con- 
dition number C for the new operator, as well 
as the decay constant 9 in the locality bound l(2Jl 
for the corresponding overlap Dirac operator, and 
compare with the values for other operators. The 
9 for the new operator is significantly larger than 
for the previously considered free field hypercubic 
operators, and is larger than the Wilson-Dirac 9 
by more than a factor of 3. 

Gauged hypercubic operators can be built up 
in various ways from free field operators. We 
now point out that if the free field operator is 
a good approximate GW solution then the same 
is generally true for the gauged operator when the 
gauge field satisfies the admissibility condition lO 
with small e. General ultra-local gauged lattice 
Dirac operators are polynomials in the parallel 
transporters T±^. Since these can be expressed 
in terms of the S'^ , Ci, defined earlier. A* A — 
A* A(S^j^,Cy). Inserting the spectral decomposi- 



tions S"^ = Y.a^ Sfj,a^Pf,a^ , = Y^ji^ Cyfi^P^fi^ 

(where P^ia^ , Qvji^ are the projections onto the 
eigenspaces with eigenvalues Sy^a^ , Cy^i^ , respec- 
tively) , and using the fact that in any gauge back- 
ground s^a^,Ci,f3^ S [—1, 1], one can make a con- 
nection with A* A free and derive bounds 

XrmnJ - Xrmnj{Cf - 1)K- - 0(e) < A* A 

< XmaxJ + A™„,/(C/ - 1)X_ + 0(e) (7) 

(the derivation will be given elsewhere) where 
the subscript "/" refers to the free field quan- 
tities. Here K_ = n",,;! | E.(-) (V-, 0.^) I where 
0ai---p4 = Piai ■ ■ ■ QiPi and the sum is restricted 
to the cr = {ai, . . . ,{3^} for which {tp^ 0^4') is neg- 
ative. Since E<t(+)(V', 0^^") + Ea(-)(V', 0^-0) = 
1 for all unit norm -0, and K?/', Oo-V') I 5; 1 ^""j 
we can expect K- to be of order 1 in typical 
gauge backgrounds. Then, if the free field op- 
erator is a good approximate GW solution (i.e. 
Cf ~ 1), the terms with K- are very small in 
((TJ. The bounds then imply that the condition 
number of the gauged operator is close to the 
free field condition number when e is small. Fi- 
nally we mention that hypercubic and other short 
range lattice Dirac operators which are reason- 
able approximate GW solutions in equilibrium 
gauge backgrounds have already been found in 
numerical work; see |10) and the ref.'s therein. 
I thank Wolfgang Bietenholz for discussions. 
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